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Abstract 

With Monte Carlo methods, we systematically investigate the short-time dynamics of domain- 
wall motion in the two-dimensional random- field Ising model with a driving field (DRFIM). We 
accurately determine the depinning transition field and critical exponents. Through two different 
definitions of the domain interface, we examine the dynamics of overhangs and islands. At the 
depinning transition, the dynamic effect of overhangs and islands reaches maximum, and this is an 
important mechanism leading the DRFIM model to a different universality class from that of the 
Edwards- Wilkinson equation with quenched disorder. 
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I. INTRODUCTION 



In the past years, much effort of physicists has been devoted to domain-wall dynamics in 
magnetic devices, nanomaterials, thin films, semiconductors, contact lines, and fluid invasion 
in porous media In particular, quenched randomness in ferroic materials, e.g. ultrathin 

'erromagnetic and ferroelectric films, fundamentally affects the response to an external field 
7H1C| . From a pragmatical point of view, understanding the controlled movement of domain 
walls plays an important role in developing new classes of potential non-volatile storage- 



class memories 
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12| . From a purely theoretical point of view, it is also essential 



for understanding non-equilibrium dynamics in disordered media [13Hl6|. For a dc (direct 
current) driving field, H, the domain-wall motion exhibits a depinning transition at zero 
temperature. The depinning field. He, separates the regimes of static pinning {H < He) 



and friction-limited viscous slide {H > He 
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3- 



At low, but non-zero, temperatures. 



;he sharp depinning transition is softened and a thermally activated creep state appears 
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22| . For an ac (alternative current) driving field, H{t) = HQexp{i2TT ft) , and at a 



nonzero temperature, the domain-wall motion exhibits different states and dynamic phase 
transitions, which can be classified in the so-called Cole-Cole diagrams 
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25|. 



Current theoretical approaches to domain-wall dynamics in ultrathin films are typi cally 



based on the Edwards- Wilkinson equation with quenched disorder (QEW) [l4, llSl, |26|, |27 |. 
This equation is a phenomenological model, and detailed microscopic structures and inter- 
actions of real materials are not concerned. Additionally, a self-inconsistence is puzzling at 



the depinning transition, especially for the roughness exponent ( |28l-l30j . To understand the 



domain-wall motion at a microscopic level, one needs lattice models based on microscopic 

3l| . The random-field Ising model with a driving field 



structures and interactions 
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(DRFIM) is an candidate, at least to capture robust features of the domain-wall motion, 
although it does not include all interactions in real materials. 

As an interface propagates, overhangs and islands may occur. Such a phenomenon is 
observed in many experiments for magnetic materials 32|, |33| . In the QEW equation, an 
interface is described by a single-valued elastic string, and overhangs and islands are not 



taken into account 
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35|. However, the anomalous roughness exponent {( ^ 1) has led 



ones to conjecture that an one-dimensional string necessarily develops overhangs and islands 
. In the DRFIM model, overhangs and islands are naturally created, thus the domain 
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wall is not single- valued and one- dimensional 



36 



. Numerical simulations of the DRFIM 



model show that characteristics of the domain-wall motion may depend on overhangs and 



islands 



39| . In the literatures, however, the dynamics of overhangs and islands is rarely 



referred, and its dynamic effect on the depinning transition is not identified. 



In 
rium 



;he past years much progress has been achieved in critical dynamics far from equilib- 
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43|. Although the spatial correlation length is still short in the beginning of the 



time evolution, the short-time dynamic scaling form is induced by the divergent correlat- 
ing time around a continuous phase transition. Based on the short-time dynamic scaling 
form, new methods for the determination of both dynamic and static critical exponents have 
been developed 
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44| . Since the measurements are carried out in the short-time regime. 



one does not suffer from critical slowing down. Recent activities include various applications 
and developments such as theoretical and numerical studies of the Josephson junction arrays 
and aging phenomena 
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49| . A kind of domain- wall roughening process at order-disorder 



phase transitions has also been revealed |50l-l52|. Very recently, the short-time behavior 



of the domain-wall relaxation around the depinning transition is noted in simulations and 
experiments Q, [isl, 531. 



In a recent article 



16| . based on the short-time dynamic approach, the domain- wall 



dynamics of a disordered magnetic system driven by a constant field H is investigated with 
the DRFIM model, in comparison with the QEW equation and experiments. The depinning 
transition at zero temperature is of second order, and its ordered parameter is the interface 
velocity. The transition field, static and dynamic exponents, and local and global roughness 
exponents are accurately determined for the DRFIM model. The results indicate that the 
DRFIM model does not belong to the universality class of the QEW equation |Ji)|], in contrast 



to the usual assumption 
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37|. 



However, it is unknown what mechanism leads the DRFIM model to a different universal- 
ity class from the QEW equation. The fluctuation and correlation of the interface velocity, 
and especially, the dynamics of overhangs and islands are not touched at all. Due to the 
existence of overhangs and islands, the definition of the domain interface is theoretically not 
unique, and comparison with experiments remains ambiguous. The purpose of this paper is 
to provide a comprehensive understanding of the depinning transition in the DRFIM model. 
We emphasize that the dynamics of overhangs and islands plays a key role. In Sec. II, the 
model and scaling analysis are described. In Sec. Ill, numerical results are presented. Sec. 
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IV includes the conclusions. 



II. MODEL AND SCALING ANALYSIS 



A. Model 



The random-field Ising model is defined by the Hamiltonian 

n = -jY. s,s,-HY,s,-Y.KSi, (i) 

<ij> i i 

where Si = ±1 is the Ising spin of the two-dimensional square lattice. The random field hi 
is uniformly distributed within an interval [—A, Al, and H is a homogeneous driving field. 

nn 

Following Refs. [16|, 117[, we fix A = 1.5 J and set J = 1. For comparison, A = is also 
simulated. A Gaussian distribution of the random field hi leads to similar results, but it is 
technically more complicated. Therefore, the results of the uniform distribution of hi are 
presented in this paper. Our simulations are performed at zero temperature with lattice sizes 
L = 128,256,512, and 1024 up to tmax = 2000, with total samples 10,000,40,000,50,000, 
and 30, 000, respectively. Simulations of different L confirm that our results do not suffer 
from finite-size effects. Errors are estimated by dividing the samples into three or four 
subgroups. If the fluctuation of the curve in the time direction is comparable with or larger 
than the statistical error, it will be taken into account. 

The initial state is such a state, that spins are positive in the sublattice on the left side 
and negative on the right side. We set the x axis in the direction perpendicular to the perfect 
domain wall. Antiperiodic and periodic boundary conditions are used in x and y directions, 
respectively. To eliminate the pinning effect irrelevant for disorder, we rotate the square 
lattice such that the initial domain wall orients in the (11) direction of the square lattice, as 



shown in Refs. 



11, 



21 



39| . After preparing the initial state, we randomly select a spin, 
and flip it if the total energy decreases after flipping. A Monte Carlo time step is defined 
by single-spin selects. As time evolves, the domain wall moves and roughens while the 
bulk remains unchanged. 

We should emphasize that overhangs and islands are naturally created in the DRFIM 
model during the time evolution, but only by the domain wall, since the temperature is set 
to zero. Therefore, it makes sense to deflne the domain wall as an interface, the so-called 
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domain interface. Due to the existence of overhangs and islands, however, there may be 
different ways to define the domain interface. In this paper, two typical definitions of the 
domain interface are concerned, to investigate the dynamics of overhangs and islands, and to 
compare with the QEW equation and experiments. We first study the definition commonly 
used in the literatures, i.e., the one define by the magnetization 
spin at site {x,y) by Sxy{t), we introduce a line magnetization 
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2l|. Denoting a 



'm{y,t) = - 



lx=l 



(2) 



The height function of the domain interface is then defined as 



hm{y,t) = -[m{y,t) + 1]. 



(3) 



The subscript m indicates that it is defined by the line magnetization. Thus the roughness 
function is introduced to depict the roughening of the domain interface, 



(4) 



where < • • • > includes the statistical average and average over y. A more informative 
quantity is the height correlation function j34|. 



Cm{r, t) = { [hm{y + r,t) - hm{y, t)f 



(5) 

It describes both the spatial correlation of the height function in the y direction and the 
growth of the domain interface in the x direction. To independently estimate the dynamic 
exponent z, we introduce an observable 



FUt) = [M'^'\t)-M{tr]/uUt). 



(6) 



Here M(t) is the (?/o6a/ magnetization and M^'^\t) is its second moment. In fact, Fm,{t) is 
nothing but the ratio of the planar susceptibility and line susceptibility. 

On the other hand, the local interface velocity is defined as the time derivative of the 

n 

height function [17J 



dhm{y,t) 
dt 



(7) 



The average velocity of the domain interface is then obtained. 



(8) 



Here we use the subscript M to emphasize that Vj^j(t) is the global average velocity. In 
fact, f^^(t) is the order parameter of the depinning phase transition. The local and global 
fluctuations of the local interface velocity, f^-'(t) and v^^\t), are interesting observables 



v^^\t) = {vUy,t)')-v,,{tr, (9) 

2, 



1 ^ 



VMitf- (10) 



In Eq. (fTOl) . < ■■■ > only includes the statistical average. We note that f^-*(t) describes 
the fluctuation of the interface velocity in the x direction, while v'^^^it) also includes the 
correlation of the interface velocity in the y direction. 

In the definition of the height function in Eq. ([3]), overhangs and islands looks formally 
suppressed. However, they do affect the dynamic evolution of the spin configuration, and 
the dynamic effect of overhangs and islands is partially included in the interface propagation 
and growth [16]. Due to the existence of overhangs and islands, however, the definition of 
the domain interface is not unique. Whether different definitions lead to the same results, 
and how to compare with the QEW equation and experiments, remain ambiguous. For 
example, another definition of the height function can be introduced by the envelop of the 
positive spins, as shown in Fig. [1] In other words, we define the height function he{y,t) as 
the largest x coordinate of the positive spins with the fixed y coordinate and at the time t. 
This definition is closer to the experiments with imaging technique [2], 32], 33 1. 

With he{y,t), one can derive the local interface velocity Ve{y,t). Similar to Eqs. ([1])- 
( ITOl) . the average velocity Vj^(t), the roughness function (jj1(t), the function Fe(t), the height 
correlation function Ce(r, if;), and the local and global fluctuations of the velocity, vf\t) and 
f^^)(t), can be calculated. In the definition of the height function with the envelop, the 
dynamic effect of overhangs and islands is maximally taken into account. In contrast to it, 
the height function defined with the magnetization includes only a minimal contribution of 
overhangs and islands. 

To reveal the dynamic characteristics of overhangs and island, therefore, we may introduce 
the overhang height function 6h{y,t) and local overhang velocity 6v{y,t) 

Sh{y, t) = Kiy, t) - h^iy, t), (11) 

6v{y,t) = Ve{y,t) - Vm{y,t). (12) 
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Then the average overhang velocity is 



(13) 



and it is related to the average overhang size by 



6h{t) 



6v{t')dt'. 



(14) 



Similar to Eqs. (jl]), and ( fTOjl . the roughness function of overhangs and islands, 5h^'^\t), 
the local and global fluctuations of the overhang velocity, 6vi (t) and 6vq (t), can be defined. 



5h^'\t) = {[sh{y,t)r)-[5hit)r, 

6vP{t) = {[6v{y,t))r)-Mt)?: 



y=l 



Mt)f 



(15) 
(16) 
(17) 



To measure the dynamic exponent independently, we can construct a function Fs{t) in a 
similar form of Eq. ([6]) . 



B. Scaling analysis 



Since the depinning transition is a second-order phase transition, the dynamic evolution 
of the order parameter Vj^^it) should obey the dynamic scaling theory supported by the 
renormalization-group calculations 

nonequilibrium correlation lengt h ^ ~ t^/^, scaling arguments lead to a dynamic scaling form 
for the order parameter 



40, m 



44j . For a finite lattice size L, and assuming a 



43 



44| 



v^it,T,L)=b-^/''Gib-%b'^^T,b-'L). 



(18) 



Here b is an arbitrary rescaling factor, /3 and u are the static exponents, z is the dynamic 
exponent, and r = (H — He) /He- Setting b ^ ^{t) ~ t^/^, the dynamic scaling form is 
rewritten as 

v,,{t,T,L) = t-^/''^G{l,t'/''^r,t-'/^L). (19) 

In the short-time regime, i.e., the regime with $^{t) ~ t^^^ ^ L, the finite-size effect is 
negligibly small, 

v^{t,T)=t-^/''^G{f/^'T). (20) 



Therefore, at the transition point r = 0, a power law behavior is obtained, 

v,,{t) = r^/^\ (21) 

With Eq. (!20|) . the critical field He may be located by searching for the best power-law 
behavior of v^j(t,T) 44 1. The critical exponent fS/uz is then estimated from Eq. ( 12T1) . 
The critical exponent 1/;/^ can be obtained from the time derivative of f^^(t,r), calculated 
according to Eq. (120|) [l6|. 

In general, even at the transition point, the roughness function uj^{t), the height corre- 
lation function Cm{r,t) and the global fiuctuation v^^\t) of the interface velocity in Eqs. 
dl]), ([5]) and f lTOj) . may not obey a perfect power-law behavior in early times. In fact, the 
domain interface and its velocity also roughen even without disorder (A = 0), due to the 
random updating scheme in numerical simulations. This may induce corrections to scaling. 
To capture the dynamic effect of disorder, we introduce the pure roughness function Duj'^{t), 
the pure height correlation function DCm{r,t), and the pure global fiuctuation Dv^^\t) by 
subtracting the contributions of A = 0, u'^^fj{t), Cm,b{r,t), and t>[^|f,(t), respectively. For a 
sufficiently large lattice and at the transition point, we should observe standard power- law 



scaling behaviors 



and 
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5l| 



(22) 



DCm{r,t) 



if r < ^{t) < L 
if < ^{t) < r 



(23) 



Here ^(t) ~ t^/^ C is the global roughness exponent, and Qoc is the local one. Finally, the 
dynamic exponent z is independently determined by 



(24) 



Compared with the velocity itself, the global and local fiuctuations of the velocity exhibit 
more complicated dynamic behaviors. In fact, there are two intrinsic length scales at the 
depinning transition, the correlation length ^ (t) ~ t^^^ in the y direction, and the character- 
istic length of roughening in the x direction, l(t) ~ t^^^, defined in Eq. ( l22l) . It only happens 
that l{t) is irrelevant for the average velocity Vf^_j{t). In general, therefore, the pure global 



fiuctuation of the velocity should obey the scaling form 43|, |44| 



Dv^^\t) = t-'^i^^Fm/LAt)/m)- 



(25) 



For a sufficiently large lattice, finite-size scaling analysis leads to 43|, |44| 



Dv^^Ht) ~ m/LY, (26) 

where ci = 1 is the spatial dimension of the domain interface. Our numerical simulations 
show that the dependence of Dv^^\t) on l(t)/^{t) takes also a power law. Therefore, 

Dv'-^^ (t) ~ t [i-^/^/'^+^i /L, (27) 

where A^^ is an exponent, reflecting the dynamic effect of the roughness of the domain 
interface. Similarly, the local fluctuation of the velocity corresponds to = 0, 

i;t2)(t)~tI-2'3/-+A™(C-i)lA. (28) 

In fact, our numerical simulations yield A^^ ~ 2 and Am ~ 3. For the standard order-disorder 
phase transition, such terms described by A^,^ and Am do not exist, since the roughness 



exponent C = 1 IS2|. This is a difference between the depinning transition and the 
standard order-disorder phase transition. 

Due to the dynamic effect of overhangs and islands, it is unclear whether the domain 
interface defined by the envelop of the positive spins also obeys the dynamic scaling forms 
described in Eqs. f|T8|) - fl28|) . In fact, we may also examine the dynamics of overhangs and 
islands independently. Similar to the determination of the critical field He from Eq. ( 120|) . one 
can locate the critical field He by searching for the best power-law behavior of the overhangs 
velocity 6v{t,H). At the critical point He, however, the overhang velocity 6v{t) increases 
with time, different from the interface velocity v{t). For example, we may write 

6v/v ~ t^. (29) 

Here 6 is estimated to be about 0.5 in our numerical simulations. In general, 6h{y,t) and 
6v{y,t) are dynamic variables independent of hm{y,t) and Vm{y,t), and one needs another 
set of critical exponents to describe their dynamic behaviors. Similar to Du'^'^^t) in Eq. fl22|) . 



the pure roughness function of overhangs and islands obeys [16|, |5l|, [52 1 



D6h^^\t) ^t^^'/'', (30) 

where (s and zs are the roughness exponent and dynamic exponent of overhangs and islands, 
respectively. The dynamic exponent zs can be independently determined from 

Fs{t) ~ t'/^'/L. (31) 



Since the overhang velocity Sv{y,t) does not show correlation in the y direction, we assume 

(2,) (2) 

that both the local and global fluctuations of the overhang velocity, 6vi (t) and 6vq {t) 
defined in Eqs. (fTB]l and (IT7|) obey 

Sv^^\t) r^f''. (32) 

On the other hand, the global fluctuation Svq"* ~ L~^, while the local one ^wp"* and the 
roughness function 5h^'^\t) are L-independent. 



III. NUMERICAL SIMULATIONS 

In Ref. y, the dynamic scaling behaviors from Eq. fl2Up up to Eq. have been 
carefully examined for the domain interface defined with the magnetization in the DRFIM 
model. The relevant critical exponents are accurately determined, and are summarized in 
Table [U in comparison with those of the QEW equation. Our first task in this paper is to 
investigate whether these dynamic scaling forms hold also for the domain interface defined 
with the envelop of the positive spins. 

In Fig. [2|^a), the average interface velocity r) is displayed for different driving field 
H. It drops rapidly down for smaller H, while approaches a constant for larger H. Due to 
the dynamic effect of overhangs and islands, f ^ {t, r) does not exhibit a power-law scaling 
behavior at the transition field He = 1.2933 located from f^^(t,r). By searching for the 
best power-law behavior of f^()f:, r), however, one could detect an alternative transition 
field He = 1.2913(4), which is slightly smaller than H^. = 1.2933(2) obtained from v^j{t,T). 
Although the difference looks small, it is not a statistical error. We believe that this difference 
is due to the dynamic effect of overhangs and islands at the non-stationary stage of the 
dynamic evolution. Anyway, the critical exponent (5/vz = 0.210(2) is obtained from the 
slope of the curve at H^ according to Eq. fl2T|) . In a similar way, we measure l/pz = 0.76(3), 



2C/z = 1.78(1) and 1/z = 0.778(7), based on Eqs. ([20]), ([22]) and ^ respectively [l6|. We 
then calculate the critical exponents (3 = 0.278(4), u = 1.02(4), z = 1.28(1) and C = 1.14(1). 

In Fig. [21(b), the pure height correlation function DCe{r,t) is displayed as a function of 
r at He for different time t . According to Eq. ([231) . the critical exponent 2(ioc = 1.13(2) is 
derived from the slope of the curve at a large time t = 2000. To fully confirm the scaling 
form of DCe{r,t) in Eq. ([23]) . for example, we fix t' = 1024, and rescale r of another t to 
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{t'jtfl^r, and DCe{r,t) to {t'/ty</'DCeir,t). As shown in Fig. [Sl^b), data of different t 
nicely collapse to the curve of t' = 1024 with ( = 1.14 and z = 1.28 as input. Hence the 
scaling form is validated. In experiments, the critical exponent doc is usually measured from 
the power-law behavior DCe{r,t) ~ r^Cioc gj^ ^ large t. Indeed, this power-law behavior of 
DCe{r,t) is much cleaner than that of DCm{r,t), as shown in the inset. Plotting DCe{r,t) 
as a function of t for different r, we measure 2(/z = 1.78(1) and 2{( — Ciocj/ z = 0.890(5), 
based on Eq. (12^ . With the dynamic exponent z = 1.28(1) as input, we may calculate 
( and (loc- The global roughness exponent ( = 1.14(1) is the same as that obtained from 
the roughness function, and the local roughness exponent C/oc = 0.569(6) is consistent with 
Cioc = 0.565(10) measured directly in Fig. [2]^b). 

All the measurements of the critical exponents and transition field are summarized in 
Table [H in comparison with those for the domain interface defined with the magnetization 
and for the QEW equation. These results further confirm that the DRFIM model and QEW 
equation are not in a same universality class 16|, |36|. For the domain interface defined with 



the magnetization, the exponents /3, z and ( of the DRFIM model differ from those of the 
QEW equation by about 10 percent, and especially, the difference of u and between two 
models reaches nearly 30 percent fl6|. For the domain interface defined with the envelop, 
the difference is even larger. For the exponent (loc, for example, the difference between two 
models is about 45 percent. These results suggest that it is mainly the overhangs and islands 
that induces the difference between the DRFIM model and QEW equation. Although real 
materials may include more complicated interactions than the DRFIM model, experimental 
measurements of the local roughness exponent of the domain interface support doc < 1- For 
T > and < H < He, for example, it is reported that Qoc = 0.7(1) and 0.69(7) in the 
experiments with ultrathin Pt/Co/Pt films [2, l7|, and C«oc = 0.78(1) with Co28Pt72 alloy 
films [33 1 . 

In the stationary state, the depinning transition should be uniquely defined. Therefore, 
the transition field located from the short-time dynamic behavior of f ^ {t, r) is an effective 
one. In other words, Vj^{t,T) is not a "good" order parameter in rigorous sense, although 



it is closer to the experiments with imaging technique [2, l32, ISSj. To demonstrate the 
difference between the domain interfaces defined with the magnetization and envelop, we 
may consider the fluctuation of the interface velocity. In Fig. [3]^a), the global fluctuations of 
the interface velocity, v^^\t) and vf^f^{t) for A = 1.5 and 0, and the pure global fluctuation 
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Dv^^\t) = v^^'^{t) — f^^\(t) are displayed. Obviously Dv^^^{t) shows a cleaner power-law 
behavior than v^^^{t) does, due to the subtraction of The slope of the curve Dv'^^it) 

is 0.532(5). Including a power-law correction, e.g., Dv^'^\t) ~ t^{l + c/t), may improve 
the fitting to the numerical data, and the resulting exponent is consistent with 0.532(5) 
within errors. According to Eq. f l271) . one can calculate the exponent A^,,^ = 2.04(5) from 
[1 - 2/3/Z/ + A,,(C -l)]/z = 0.532. In Fig. EJ^b), the local fluctuation v^^}{t) is plotted, and 
the slope of the curve is 0.112(3). According to Eq. ( |28|) . A^,^ = 3.06(3) can be derived from 
[-2/3/z/ + A^(C - l)]/z = -0.112. To study possible finite-size effects, v^^\t) and t;^)(t) 
computed with different lattice sizes at He = 1.2933 are also shown in Fig. El All the curves 
in Fig. [3](a) are rescaled by a factor L, based on Eq. ( 1271) . As it can be seen in the figure, 
the finite-size effect can be easily controlled, i-e., it drops rapidly as L increases. This is a 
merit of the short-time dynamic approach |l6l. 1421444 1. 

However, both Dv^^\t) and v^\t) do not exhibit a power-law behavior either at He = 
1.2913 or 1.2933. In other words, they do not obey the scaling forms in Eqs. (1271) and 
( l28l) . Additionally, Ve{y,t) does not show a standard correlation in the y direction. In 
fact, the domain interface defined with the envelop may be considered as adding overhangs 
and islands to the domain interface defined with the magnetization. To understand the 
domain interface defined with the envelop, we may alternatively investigate the dynamics of 
overhangs and islands. In the inset of Fig. Hl^a), we present the overhang velocity 6v(t) for 
different driving field H. Different from the velocity fj,^(t) and Vj^{t) of the domain interface, 
6v{t) drops rapidly down for both smaller and larger driving field H, and reaches maximum 
at the transition point He, which is estimated to be between 1.29 and 1.30. In order to 
determine the transition field more accurately, we calculate the overhang size 6h{t, H) as 
a function of the external field H dX t = 1000. As shown in Fig. 111(a), the maximum of 
Sh{H) yields He = 1.294(1), in good agreement with the transition field He = 1.2933(2) 
estimated from Vj^j{t). The dynamic effect of overhangs and islands is the most prominent 
at the depinning transition point. In Fig. 111(b), 5v/v is plotted a.t H = 1.2933 with different 
lattice size L. The finite-size effect is negligible for L = 1024 up to t = 2000. According to 
Eq. ( l29l) . 6 = 0.50(2) is obtained with a power-law fit for t > 100, and a power-law correction 
improves the fitting to the numerical data. 

To determine the dynamic exponent zs, we plot the function Fs{t) at He = 1.2933 in 
Fig. I5]^a). The slope of the curve with L = 512 is 0.90(1), and it yields zs = 1.11(1). An 
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power-law correction may improve the fitting to the numerical data. Collapse of the curves 
with L = 256 and 512 indicates the L-dependence Fs{t) ~ 1/L. In Fig. Mjo), the roughness 
functions 5h'^'^\t) and 6h^^\t) for A = 1.5 and 0, and the pure roughness function D6h^'^\t) 
are displayed at H = 1.2933. Obviously, D5h^'^\t) exhibits a cleaner power-law behavior 
than Sh^'^\t). The slop of the curve is 2.10(3), and it leads to Cs = 1.16(2), close to C = 1.14 
for the domain interface. With a correction to scaling, i.e., D5h^'^\t) ~ t'^'^^^^^{l + c/t"^), the 
fitting to numerical data is extended. Since both the overhang size and velocity reaches a 
maximum at the transition field He, the standard exponent u is not defined. In other words, 
dnSvit, H) is close to zero at H^, and effectively, vs ^ 1- 

Finally, we consider the local and global fluctuations of the overhang velocity. In Fig. 0, 
the local fluctuation 5vf'\t) is plotted at H = 1.2933. By the definition, the local fluctuation 
can also be calculated with 6vi'^\t) = vf\t)—v''^{t) — 2/\{t), and A(t) =< Vm{y,t)6v{y,t) > 
—Vj^,j{t)Sv{t). Hence v^^\t) — vl^\t) is plotted for comparison. The result indicates that 
6vl'^\t) ~ v^^\t) — v^{t), i.e., /S.{t) ~ 0. According to Eq. ( l32l) . a direct measurement 
from the slope of the curve gives 2a = 0.972(4). A power-law correction to scaling yields 
a similar result 2a = 0.976. The global fluctuation 6vQ\t) and v^^\t) — v^^\t) are also 
plotted in Fig. [61 Both of them are rescaled by a factor of L/A, because of the finite-size 
dependence of L~^. Overlapping of these two curves is observed and it yields A(t) ^ 0, 
too. The exponent 2a = 1.002(5) is determined from the slope, in agreement with that 
measured from the local fluctuation of the velocity. Actually, our numerical simulations 
show that 5vi (t) ~ L Svq (t). In Fig. [6], it is only for clarity that 5vq (t) is rescaled by 
a factor of L/4 rather than L. The result A(t) ^ indicates that Vm{y,t) and 6v{y,t) are 
not correlated, and therefore, another set of critical exponents is needed for describing the 
dynamic behavior of overhangs and islands. 

The fact, that both 6vi (t) and 6vq (t) are governed by a same exponent a ~ 0.5, 
indicates that the overhang velocity 6v{y,t) is not correlated in the y direction, although 
the overhang size 6h{y,t) does. If we consider 6v{y,t) as "a height function", its roughening 
process described by 6vi (t) belongs to the universality class of random depositions. Why 
does such a phenomenon occur? For example, he{y + l,t) may suddenly produce an overhang 
at he{y,t), and induce a rapid increase of the overhang velocity 6v{y,t). However, this is 
not correlated with 6v{y + l,t). For the domain interface defined with the envelop, Ve{y,t) 
is also not correlated in the y direction, since Ve{y,t) = Vm{y,t) + 6v{y,t), and 6v{y,t) 
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dominates at larger t. At the transition field He, our numerical simulations show that 
vP{t) ^ vl^^t) + Svl^\t), and v'^^\t) ^ v'^\t) + dv^^^t). On the other hand, Svi^\t) and 
^'^c ^ (t) exhibit different power-law behaviors from (t) and v^^^ (t) respectively. Therefore, 
f (t) and f (t) do not obey a simple power law. For large t, however, 5f (t) and 5f (t) 
dominate the dynamic behaviors of v^^\t) and v^^\t). 



IV. SUMMARY 



Based on the short-time dynamic approach, we have systematically investigated the 
domain-wall dynamics of the DRFIM model at the depinning transition, and have accu- 
rately determined the transition field and all the static and dynamic critical exponents. 
Through two different definitions of the domain interface, we examine the dynamics of over- 
hangs and islands. All the critical exponents for the domain interface and for overhangs and 
islands are summarized in Table [H 

For the domain interface defined with the envelop, we do observe the dynamic scaling 
behaviors in Eqs. (!20|) - fl2^ at an effective transition field He = 1.2913, slightly below He = 
1.2933 determined from the interface velocity defined with the magnetization. The difference 
of the critical exponents between the DRFIM model and QEW equation becomes lager for 
the domain interface defined with the envelop, especially for the local roughness exponent 
Cioc- These results further support that the DRFIM model and QEW equation are not in 
a same universality class. Since the dynamic effect of overhangs and islands is maximally 
expressed in the domain interface defined with the envelop. We argue that it is mainly the 
overhangs and islands that induces the difference between the DR! 
equation. Experiments report the local roughness exponent Qoc < 1 

The global and local fluctuations of the interface velocity deflned with the magnetization 
exhibit the power-law scaling behaviors described by the exponents Xj^^ and Am in Eqs. fl27]) 
and fl28|) . while those deflned with the envelop do not. It indicates that the interface velocity 
deflned with the envelop is not a good order parameter of the depinning transition, although 
it is closer to experiments with imaging techniques. In fact, Ve{y,t) is also not correlated in 
the y direction. It should be interesting to measure the fluctuation and correlation of the 
interface velocity in experiments. 

At the depinning transition, the dynamics effect of overhangs and islands reaches max- 
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imum. The observables of overhangs and islands do obey the dynamic scahng forms in 
Eqs. f l29|) -f l32|) . similar to those for the domain interface. However, another set of critical ex- 
ponents should be introduced, since Vm{y,t) and 6v{y,t) are not correlated. Different from 
the interface velocity, the overhang velocity increases with time. The dynamic exponent 
zs = 1.11(1) for overhangs and islands is smaller than zs = 1.33(1) for the domain interface, 
while the roughness exponent is about the same. In particular, the overhang velocity 6v{y, t) 
does not show correlation in the y direction, although the overhang size 5h{y, t) does. Con- 
sidering Sviy, t) as a height function, its roughening process belongs to the universality class 
of random depositions. Since hf,{ii, t) = hm{y, t) +6h{y, t), the domain interface defined with 
the envelop can be mostly understood by adding overhangs and islands to that defined with 
the magnetization. 
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QEW 


Magnetization 


Envelop 


Overhang 


vit) 


He 




1.2933(2) 


1.2913(4) 


1.294(1) 




fi 


0.33(2); 0.33 


0.295(3) 


0.278(4) 






V 


1.29(5); 1.33; 1.33(1) 


1.02(2) 


1.02(4) 


» 1 




z 
9 


1.5; 1.53 


1.33(1) 


1.28(1) 


1.11(1) 
0.50(2) 




c 


1.26(1); 1.25; 1.24 


1.14(1) 


1.14(1) 


1.16(2) 


Cir,t) 


c 


1.23(1); 1.25 


1.13(1) 


1.14(1) 






Cloc 


0.98; 0.92 


0.735(8) 


0.569(6) 




/'o^ / \ 


A 




2.04(5) 












3.06(3) 








a 








0.501(3) 












0.488(4) 



TABLE I: The depinning transition field and critical exponents obtained for the DRFIM model 

3, Q: 



are compared with those for the QEW equation in Refs. 
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54l |. The exponerits ^, v, 



z, C and Ooc for the domain interface defined with the magnetization are taken from Ref. [16iJ 
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(a) (b) 

FIG. 2: (a) Interface velocity v^{t,T) is plotted with solid lines for different driving fields H with 
L = 512 on a log-log scale. For clarity, the curve of H = 1.2913 is shifted down. For comparison, 
the curve with L = 1024 is shown with open circles, (b) The pure height correlation function 
DCe{r,t) is displayed at He = 1.2913. According to Eq. (p6]) . data collapse is demonstrated. 
Stars, triangles, pluses, and circles correspond to t = 128, 256, 512, and 2000, respectively. In the 
inset, DC{r,t) defined with the envelop (upper) and with the magnetization (lower) are shown at 
t = 2000. In both (a) and (b), dashed lines show power-law fits. 
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1000 



FIG. 3: In (a) and (b), the global and local fluctuations of the interface velocity are displayed at 
He = 1.2933 for different L on a log-log scale. Dashed lines show power-law fits, and the solid line 
represents a power-law fit with correction. In (a), the pure global fluctuation function Dv^^\t) is 
shown for L = 1024 with squares. All the curves have been rescaled by a factor L. 




(a) (b) 

FIG. 4: (a) The overhang size 6h{t, H) is plotted as a function of the driving field H at t = 1000. 
In the inset, the overhang velocity 5v{t) is shown for different i? on a log-log scale. For clarity, the 
curves for H = 1.29, 1.28 and 1.27 are shifted up. (b) 6v/v is displayed at = 1.2933 for different 
L on a log-log scale. The solid line represents a power-law fit with correction. 
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(a) (b) 

FIG. 5: (a) Fs{t) is plotted at He = 1.2933. The curve of L = 256 has been rescaled by a factor 
1/2, according to Eq. (j3ip . (b) The roughness function and pure roughness function of overhangs 
and islands are displayed at He = 1.2933. In both (a) and (b), dashed lines show power-law fits, 
while solid lines are for power-law fits with correction. 




10 100 1000 

t 



FIG. 6: The local and global fluctuations of the overhang velocity, dvl and Svq , are plotted on a 
log-log scale. For comparison, — and v^^^ — v^^^ are also displayed. To show the finite-size 
dependence, the global quantities are rescaled by a factor of L/4. Dashed lines show power-law 
fits, while solid lines are for power-law fits with correction. 
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